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Abstract 

We present a rigorous derivation of the flow at arbitrary time in a deterministic 
cellular automaton model of traffic flow. The derivation employs regularities in preim- 
ages of blocks of zeros, reducing the problem of preimage enumeration to a well known 
lattice path counting problem. Assuming infinite lattice size and random initial con- 
figuration, the flow can be expressed in terms of generalized hypergeometric function. 
We show that the steady state limit agrees with previously published results. 

1. Introduction 

Since the introduction of the Nagel-Schreckenberg (N-S) model in 1992 [8], cellular automata 
became a well established method of traffic flow modeling. Comparatively low computational 
cost of cellular automata models made it possible to conduct large-scale real-time simulations 
of urban traffic in the city of Duisburg [2] and Dallas/Forth Worth [10]. Several simplified 
models have been proposed, including models based on deterministic cellular automata. For 
example, Nagel and Herrmann [9] considered deterministic version of the N-S model, while 
Fukui and Ishibashi [5] introduced another model (to be referred to as F-I model), which 
can be understood as a generalization of cellular automaton rule 184. Rule 184, one of the 
elementary CA rules investigated by Wolfram [13], had been later studied in detail as a 
simple model of surface growth [6], as well as in the context of density classification problem 
[3]. It is one of the only two (symmetric) non-trivial elementary rules conserving the number 
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of active sites [1], and, therefore, can be interpreted as a rule governing dynamics of particles 
(cars). Particles (cars) move to the left if their right neighbor site is empty, and do not move 
if the right neighbor site is occupied, all of them moving simultaneously at each discrete 
time step. Using terminology of lattice stochastic processes, rule 184 can be viewed as a 
discrete-time version of totally asymmetric simple exclusion process. Further generalization 
of the F-I model has been proposed in [4]. 

In all traffic models, the main quantity of interest is the average velocity of cars, or the 
average flow, defined as a product of the average velocity and the density of cars. The graph 
of the flow as a function of density is called a fundamental diagram, and is typically studied 
in the steady state it — > oo). For the F-I model, steady-state fundamental diagram can be 
obtained using mean-field argument [5], as well as by statistical mechanical approach [11] or 
by studying the time evolution of inter-car spacing [12]. In general, little is known about non- 
equilibrium properties of the flow. In [3], we investigated dynamics of rule 184 and derived 
expression for the flow at arbitrary time, assuming that the initial configuration (at t = 0) 
was random, using the concept of defects and analyzing the dynamics of their collisions. In 
what follows, we shall generalize results of [3] for the deterministic F-I traffic flow model and 
derive the expression for the flow at arbitrary time. The derivation employs regularities of 
preimages of blocks of zeros, reducing the problem of preimage enumeration to a well known 
combinatorial problem of lattice path counting. Assuming infinite lattice size and random 
initial configuration, the flow can then be expressed in terms of generalized hypergeometric 
function. We will, unlike in [3], explore regularities of preimages using purely algebraic 
methods, i.e., without resorting to properties of spatiotemporal diagrams and dynamics of 
defects. 

2. Deterministic traffic rules 

Deterministic version of the F-I traffic model is defined on one-dimensional lattice of L sites 
with periodic boundary conditions. Each site is either occupied by a vehicle, or empty. The 
velocity of each vehicle is an integer between and m. If x(i,t) denotes the position of 
the ith car at time t, the position of the next car ahead at time t is x(i + l,t). With this 
notation, the system evolves according to a synchronous rule given by 



is the velocity of car i at time t. Since g = x(i + l,t) — x(i,t) — 1 is the gap (number of 
empty sites) between cars i and i + 1 at time t, one could say that each time step, each car 
advances by g sites to the right if g < m, and by m sites if g > m. When m = 1, this model 
is equivalent to elementary cellular automaton rule 184, for which a number of exact results 
is known [6, 3]. 

The main quantities of interest in this paper will be the average velocity of cars at time 
t defined as 





where 






(3) 
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and the average flow (f)(t) = pv(t), where p = N/ L is the density of cars. In what follows, we 
will assume that at t = the cars are randomly distributed on the lattice. When N — > oo, 
this corresponds to a situation when sites are occupied by a car with probability /?, or are 
empty with probability 1 — p. 

In general, if Nk(t) is the number of cars with velocity k, we have 

-i m 

W) = ^T.kN k (t). (4) 

k=l 

When k < m, Nk(t) is just the number of blocks of type 10 fc l, where k denotes k zeros. 
This means that a probability of an occurrence of the block 10 fc l at time t can be written 
as P t (10 fc l) = Nk/L. Similarly, for k = m, P t (10 m ) = N m (t)/L. As a consequence, equation 
(4) becomes 

-H.p,(10M) +m P,(10-) 

We will now demonstrate that in the deterministic F-I model with maximum speed m the 
average flow depends only on one block probability. More precisely, we shall prove the 
following: 

Proposition 1. In the deterministic F-I model with the maximum speed m, the average 
flow 4> m (t) is given by 

<j) m (t) = l-p-P t (O m+1 ). (6) 



To prove this proposition by induction, we first note that for m = 1 equation (5) gives 
4>i(t) = Pf(10). Using consistency condition for block probabilities P*(10) + Pt (00) = -Pt(O) = 
1 — /?, we obtain <f>i(t) = 1 — p — P(00), which verifies (6) in the m = 1 case. Now assume 
that (6) is true for some m = n — 1 (where n > 1), and compute 4> n (t): 

n-l 

Mt) = nP t {io n ) + J2jP{^n) = 

n-2 

= nP t {10 n ) + {n- l^lO"" 1 !) + ]T jP{Wl) 

n-2 

= (n - l^lO"- 1 !) + P t (10 n )] + P t (10 n ) + i^(10 j l) 

3 = 1 

Using consistency condition P t (10 n_1 l) + P t (10 n ) = P t (10 n_1 ) we obtain 

n-2 

Taking into account that P t (10 n ) = Pt(0 n ) — P t (0 n+1 ) (which, again, is just a consistency 
condition for block probabilities), and using (6) to express </> m _i(t), we finally obtain 

</> m (t) = l- /O -P t (0 m+1 ). (7) 

This means that validity of (6) for m = n follows from its validity for m = n — 1, concluding 
our proof by induction. 
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3. Enumeration of preimages of 



Proposition 1 reduces the problem of computing (f> m (t) to the problem of finding the probabil- 
ity of a block of ra + 1 zeros. In order to find this probability, we will now use the fact that the 
deterministic F-I model is equivalent to a cellular automaton defined as follows. Let s(i,t) 
denotes the state of a lattice site i at time t (note that i now labels consecutive lattice sites, 
not consecutive cars), where s(i,t) = 1 for a site occupied by a car and s(i,t) = otherwise. 
We can immediately realize that if a site i is empty at time t, then at time t + 1 it can become 
occupied by a car arriving from the left, but not from a site further than i — m. Similarly, 
if a site i is occupied, it will become empty at the next time step only and only if site i + 1 
is empty. Thus, in general, s(i, t + 1) depends on s(i — m,t), s(i — m + 1, t), . . . , s(i + 1, t), 
i.e., on the state of m sites to the left, one site to the right, and itself, but not on any other 
site, what can be expressed as 



where f m is called a local function of the cellular automaton. For the F-I CA, one can write 
explicit formula 1 for / m , such as 



which, using terminology of cellular automata theory, represents a rule with left radius m 
and right radius 1. In general, after t iteration of this cellular automaton rule, state of a site 
s(i, t) depends on s(i — mt, 0), s(i — mt + 1,0),..., s(i + t, 0), but not on any other sites in 
the initial configuration. Similarly, a block of k sites s(i, t)s(i + 1, t) . . . s(i + k) depends only 
on a block s(i — mt, 0), s(i — mt + 1,0),..., s(i + k + 1, 0), as schematically shown in Figure 1. 
We will say that s(i — mt, 0), s(i — mt + 1,0),..., s(i + k + t, 0) is an t-step preimage of the 
block s(i, t)s(i + 1, t) . . . s(i + k). Preimages in the F-I cellular automaton have the following 
property: 

Proposition 2. Block aia 2 a 3 . . . a p is an n-step preimage of a block m+1 if and only if 
p = (n + l)(m + 1) and, for every k (1 < k < p) 



where £(1) = —m and £(0) = 1. 

Before we present a proof of this proposition, note that it can be interpreted as follows. 
Let us assume that we have a block of zeros and ones of length p, where p = (n + l)(m + 1), 
and we want to check if this block is an n-step preimage of a block m+1 . We start with a 
"capital" equal to zero. Now we move from the leftmost site to the right, and every time we 
encounter 0, we increase our capital by m. Every time we encounter 1, our capital decreases 

1 Since formula (9) will not be used in subsequent calculations, we give it withot proof (which is elemen- 
tary). 






k 




(10) 
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Figure 1: Fragment of a spatio-temporal diagram for the F-I rule with m = 2. States of nine 
sites during three consecutive time steps are are shown, black circles representing occupied 
sites. Block B\ = 101110100 is a 2-step preimage of the block B 2 = 100. Outlined sites 
constitute "light cone" of the block B 2} meaning that the state of sites belonging to B 2 can 
depend only on sites inside the outlined region, but not on sites outside this region. 

by 1. If we can move from ai to a p and our capital stays always larger than zero, the string 
aia 2 a 3 . . . a p is a preimage of m+1 . Condition (10) can be also written as 

X><^-T, (11) 
, m + 1 

8 = 1 1 

because £(x) = 1 — (ra + l)x for x £ {0, 1}. 

For the purpose of the proof, strings a\a 2 . . . a p of length p satisfying (11) for a given m 
and for every k < N will be called m- admissible strings. 
Lemma. Let s(l, t)s(2, t) . . . s(p,t) be an m- admissible string. If 

s(i,t + 1) = f m (s(i - m,t),s(i - m + l,t), . . . , s(i + l,t)J, (12) 

and if f m is a local function of the deterministic F-I model with maximum speed m, then 
s(m + 1, t + l)s(m + 2, t + 1) . . . s(p — 1, t + 1) is also an m-admissible string. 

To prove the lemma, it is helpful to employ the fact that the FT rule conserves the 
number of cars. Let < k < p and let us consider strings S\ = s(l,t)s(2,t) . . .s(k } t) and 
S2 = s(m + 1, t + l)s(2,t + l)...s(A;,t + l). If the string s(l, t)s(2\ t) . . . s(k y t) is m-admissible, 
then its hrst m + 1 sites must be zeros. This means that in one time step, no car can enter 
string s(l, t)s(2, t) . . . s{k, t) from the left. On the other hand, in a single time step, only one 
car (or none) can leave the string on the right hand side, i.e., 

k k 

£>(i,*) = e+ E s{i,t + l) (13) 

i=l i=m-\-l 

where e £ {0, 1}. Three cases can be distinguished: 
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(i) All sites s(k — m + 1, t)s{k — m + 2, t) . . . s(k, t) are empty (equal to 0). Then no car 
leaves Si, which means that e = 0, and 

k k k — m 7 

£>(i,f)= £ s(i,t + i)= £ 5 ^' t + 1 ) < ^TT- ( 14 ) 

8 = 1 i=m-\-l i=m-\-l 

The last inequality is a direct consequence of m-admissibility of S\. Since the length of the 
string S2 is equal to k — m, the above relation (which holds for arbitrary k) proves that S2 
is also m-admissible in the case considered. 

(ii) Among sites s(k — m + l,t)s(k — m + 2,t) . . .s(k,t) there is at least one which is 
occupied (equal to 1), and s(k + l,t) = 1. In this case, since the last site in S\ is "blocked" 
by the car at s(k + l,t), again no car can leave string S\ in one time step. Therefore, 

k k 

£>(M)= £ + 1). (15) 

i=l i=m-\-l 



m-admissibility of S2 implies 



k+ \> k f,s(i,t) = j:s(i,t)+i. (i6) 



m + ^ ■ , ■ , 

1 8 = 1 8 = 1 

Combining (15) with (16) we obtain 

E, \ kj , , 

s{ht + l )<——, 17 
, 1 m + 1 

8=m + l 1 

which again shows that S2 is m-admissible. 

(iii) Among sites s(k — m + l,t)s(A; — m + 2,t) . . .s(k,t) there is at least one which is 
occupied (equal to 1), and s(k + l,t) = 0. In this case, one car will leave right end of the 
string Si, therefore 

k k 

£>(M)= £ s(M + i)-i. (18) 

z = l «'=m + l 

As before, from m-admissibility of Si we have 

X>(M) = !>(*, t)<A±I (19) 

■ , r— ' m + 1 

8=1 8=1 1 

hence 

E, , x — ^ y k ~\~ 1 k / ? ? , , 

*»,f + l = £ s(t,t + l - l<— — -l = — — , 20 
, 1 • 1 m + 1 m + 1 

8=m + l 8=m + l 

which demonstrates that case (iii) also leads to m-admissibility of S2, concluding the proof 
of our lemma. 

Let us now assume that the block B\ = s(l,t)s(2,t) . . .s(p,t) is m admissible (n being 
some fixed integer and p = (m + l)(n + 1)). Applying the lemma to this block we conclude 
that B 2 = s(m + 1, t + l)s(2, t + 1) . . . s(p — 1, t + 1) is m-admissible as well. Applying the 
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Figure 2: m-admissible block with n zeros and n\ ones is equivalent to a lattice path from 
the origin to (n ,ni) which does not touch nor cross the line x = my. corresponds to a 
horizontal segment, while 1 to a vertical segment. 

lemma to B 2 we obtain m-admissible block B 3 = s(2m + 1, t + 2)s(2, t + 2) . . . s(p — 2, t + 2). 
After n applications of the lemma we end up with the conclusion that the string B n+ \ = 
s(nm + 1, n + l)s(nm + 2, n + 1) . . . — n) is m-admissible. Since the length of B n+ \ is 
p — n — nm = (n + l)(m + 1) — n(m + l) = m + l,it must, to be m-admissible, be composed of 
all zeros, i.e., B m+ i = m+1 . This means that m-admissibility of B\ is a sufficient condition 
for B\ to be an n-step preimage of m+1 . Reversing steps in the above reasoning, one can 
show that is is also a necessary condition. 

4. Fundamental diagram 

We shall now use proposition 2 to calculate P t (0 m+1 ). First of all, we note that P t (0 m+1 ) 
is equal to the probability of occurrence of t-step preimage of m+1 in the initial (random) 
configuration, that is 

^(O m+1 ) = E P o(a), (21) 

where the sum goes over all t-step preimages of m+1 . Consider now a string which contains 
n zeros and n\ ones. The number of such strings can be immediately obtained if we realize 
that it is equal to the number of lattice paths from the origin to (n , n\) which do not touch 
nor cross the line x = my } as shown in Figure 2. This is a well known combinatorial problem 
[7] , and the number of aforementioned paths equals 

n — mrii (riQ + n\ 
riQ + n x I m 




7 




Figure 3: Graph of the probability P t (0 m+1 ) as a function of p for m = 2 and t = 1 (upper 
line), t = 5 (middle line), and t = 100 (lower line). 



Probability of occurrence of such a block in a random configurations is, therefore, 

n — rani / + ^i 



n + rai 



ra 



(23) 



where /? = P(l). In a t-step preimage of m+1 the minimum number of zeros is 1 + mit + 1) 
zeros, while the maximum is (ra + l)(t + 1) (corresponding to all zeros). Therefore, summing 
over all possible number of zeros i } we obtain 



Pt(0 



ra + 1 ^ 



(m+W+i) ^_ m [( m + ^ + 1) - i] ^ ( m + l)(t + 1) 

i=l+m(t+l) 



( m + l)( t + l)_i 



(m + l)(t + 1) 
Changing summation index j = i — mit + 1) we obtain 



PJQ m+1 ) = V 3 A m + ^ + ^ | n t+1 ~ 3 (l - o) m(t+l)+3 



(24) 



Figure 3. shows a graph of P t (0 m+1 ) as a function of p for m = 2 and several values of 
t. We can observe that as t increases, the graph becomes "sharper" at p = 1/3, eventually 
developing singularity (discontinuity in the hrst derivative) at p = 1/3. More precisely, one 
can show (see appendix) that 



lim P t (0 

t-too 



m+1 = f 1 - (m + l)p if p < l/(m + 1), 
I otherwise. 



(25) 
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-Poo(O m+1 ), therefore, can be viewed as the order parameter in a phase transition with critical 
point at p = l/(m + 1). Using Proposition 1 we can now find the average flow in the steady 
state 

6 (oo) = i mf> if P< 1 /( m + 1 ) 5 (26) 
m 1 1 — p otherwise, 

which agrees with mean-field type calculations reported in [5] as well as with results of 
[11, 12]- 

To verify validity of the result for t < oo, we performed computer simulations using a 
lattice of 10 5 sites with periodic boundary conditions. The average flow has been recorded 
after each iteration up to t = 100 for three values of p: at the critical point p = 1/3 as well 
as below and above the critical point. The resulting plots of the flow as a function of time 
are presented in Figure 4. Again, the agreement with theoretical curves 

Mt) = i - p - E jjj ( ( Y+w " ») m(t+1)+ ' W 

is very good. 

Without going into details, we note that the formula (27) can be also expressed in terms 
of generalized hypergeometric function 2 Fi: 



, , _ _ (1 - /o) 1+m+m y(l + m + f + mf)! 

0m(j " 9 (1 + m + mt)(l + t)\{m + mt)\ 21 



2 + m + mt ' p 



(28) 



Since fast numerical algorithms for computing 2 Fi exist, this form might be useful for the 
purpose of numerical evaluation of (f) m (t). 



5. Conclusion 

We presented derivation of the flow at arbitrary time in the deterministic F-I cellular au- 
tomaton model of traffic flow. First, we showed that the flow can be expressed by the 
probability of occurrence of the block of m + 1 zeros P(0 m+1 ). By employing regularities in 
preimages of blocks of zeros, we reduced the problem of preimage enumeration to the lattice 
path counting problem. Finally, we used the number of preimages to find P(0 m+1 ), which 
determines the flow. 

We also found that the flow in the steady state, obtained by taking t — > oo limit, agrees 
with previously reported mean-field type calculations, meaning that in the case of the F-I 
model mean-field approximation gives exact results. This seems to be true not only for the F-I 
model, but also for many other CA rules conserving the number of active sites ("conservative" 
CA). For example, in [1] we reported that the third order local structure approximation, 
which is a generalization of simple mean-field theory incorporating short-range correlations, 
yields the fundamental diagram for rule 60200 (one of the 4-input "conservative" CA rules) 
in extremely good agreement with computer simulations. Taking this into account, we 
conjecture that the local structure approximation gives exact fundamental diagram for almost 
all "conservative" rules, excluding, perhaps, those rules for which the fundamental diagram 
is not sufficiently "regular" (meaning not piecewise linear). This problem is currently under 
investigation. 
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Figure 4: Plots of c/> 2 (t)/c/>2(oo) as a function of time for p = 0.3, p = 1/3, and p = 0.35 
obtained from computer simulation on a lattice of 10 5 sites. Continuous line corresponds to 
the theoretical result obtained using eq. (28). 
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Appendix 

In order to find the limit lim^oo P t (0 m+1 ) we can write eq. (24) in the form 



is the distribution function of the binomial distribution. Using de Moivre-Laplace limit 
theorem, binomial distribution for large n can be approximated by the normal distribution 




(29) 



where 




(30) 





1 



exp 




(31) 
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To simplify notation, let us define T = t-\- 1 and M = m + 1. Now, using (31) to approximate 
b(T — j, MT } p) in (29), and approximating sum by an integral, we obtain 

, , u rT x 1 — (T - x - MTp) 2 , , , 



'i TJ 2 irMTp(l-p) 2MTp(l-p) 

Integration yields 



-(1 — Mp) I erf f M " r - erf ( l ~ T + M ' r ^ 1 



2 I \y/2Mp(l-p)T) \j2Mp(l-p)T, 

where erf(x) denotes the error function 



erf(x) = 4= I" e~ t2 dt. (33) 

\/7T JO 



The first term in the above equation (involving two exponentials) tends to with T — > oo. 
Moreover, since lim^oo erf (a?) = 1, we obtain 

lim P t (0 m+1 ) = -(1 -Mp)\l- lim erf I 1 ~ T + M P T = ) I _ 



Now, noting that 



, 1-T + MpT j _ I 1, ifM„>l. 

t^ to \ . / 2 Mp(l - p)T J I _1 ' otherwise, 



and returning to the original notation, we recover eq. (25): 

lim P t (0^) = { I ~ {m + VP 1 {P < 1/{m + V> (35) 
t->oo otherwise. 
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